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ABSTRACT 


We measure the hadronic branching ratios of the boson into heavy quarks: 
Rb = r^o and Rc r^o— •tcc/'^—^hadrons usiug a uiulti tag 

technique. The measurement was performed using about 400,000 hadronic 
events recorded in the SLD experiment at SLAG between 1996 and 1998. 
The small and stable SLC beam spot and the CCD-based vertex detector 
were used to reconstruct bottom and charm hadron decay vertices with high 
efficiency and purity, which enables us to measure most efficiencies from data. 
We obtain, 

Rb = 0.21604 ± 0.00098(stot.) ± 0.00073(s|/st.) ^ 0.00012(Rc) 

and, 

Rc = 0.1744 ± 0.0031(stat.) ± 0.0020{syst.) T 0.0006(Rfe) 


* Work supported by Department of Energy contract DE-AC03-76SF00515 (SLAG). 
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1 Introduction 


The dominant production of decays, with large statistics, at exper¬ 
iments operating on the peak, provides a unique opportunity for probing 
electroweak interactions at high precision. The democratic production of 
all fermion flavors in decays with a clean initial state allows particu¬ 
larly sensitive tests of the Standard Model (SM) through measurements of 
fib = r zo^bb/^zo^hadrons and fi^ = Tzo^cc/'^zo^hadrons, the heavy quark pro¬ 
duction fractions in hadronic Z^ decays. The b and c quarks are the heaviest 
charge 1/3 and charge 2/3 quarks, respectively, that are accessible at the Z^ 
energy. The fif, measurement has been traditionally a focus of attention as it 
is not only sensitive to the heavy top quark mass, but also widely regarded 
as a promising window for detecting new physics through radiative correc¬ 
tions to the Z^ bb coupling. More generally, any unexpected difference in 
quark coupling of one flavor compared with other flavors could be a vital clue 
toward a solution to the puzzle of fermion family degeneracy. 

To take full advantage of the large decay samples as an effective flavor 
physics arena, the advance in vertex detector technologies has played a key 
role for flavor identihcation. Unprecedented performance in b quark tagging 
has made it possible to achieve Rb measurements mia at below 1% precision. 
Our previous Rb measurement [2j has demonstrated the effectiveness of the 
combination of high resolution vertexing and the small and stable interaction 
point at the SLC, which pointed the way for very high purity 6-tags and led 
to the reduction of systematic uncertainties for Rb measurements in general. 
In this paper we present an updated Rb measurement from SLD using ~2.5 
times more statistics than our previous publication and with an upgraded 
CCD pixel vertex detector, which improves the 6-tag efficiency to approxi¬ 
mately a factor of two higher than 6-tags used in existing measurements at 
similar 6-purity. 

The existing measurements of the charm branching ratio fic |3] are con¬ 
siderably less precise than the Rb measurements. These determinations of fic 
rely on a number of methods to identify charm jets that all have certain dis¬ 
advantages. Exclusive reconstruction of charmed mesons suffers from a small 
branching fraction and introduces a dependence on the actual value of this 
branching fraction and on the production fractions of the charm hadrons. In¬ 
clusive reconstruction using leptons or slow pions suffer from low purities and 
again a dependence on branching fractions. Due to the superb performance 
of the SLD vertex detector we are able to introduce an inclusive charm tag. 
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based on lifetime information in a similar way as the b tag, which combines 
high efficiency with good purity. Moreover, the tagging efficiencies are mea¬ 
sured from data and result in a minimal reliance on Monte Carlo simulation 
and a small systematic uncertainty. This leads to a new determination of Rc 
with much improved overall precision than the existing measurements. 

2 Apparatus and Hadronic Event Selection 

This analysis is based on approximately 400,000 hadronic events produced 
in e^e“ annihilations at a mean center-of-mass energy of ^/s = 91.28 GeV at 
the SLAG Linear Gollider (SLG), and recorded in the SLG Large Detector 
(SLD) between 1996 and 1998. A general description of the SLD can be 
found elsewhere jl]. The trigger and initial selection criteria for hadronic 
decays are described in Ref. u The Gentral Drift Ghamber (GDG) [H] and 
the upgraded Vertex Detector (VXDS) [Zj, inside a uniform axial magnetic 
field of 0.6T, provide the momentum measurements of charged tracks and 
precision vertex information near the interaction point (IP), which are central 
to this analysis. The energies of clusters measured in the Liquid Argon 
Galorimeter [S] are used for event selection and calculation of the event thrust 
axis. 

SLD uses a coordinate system with the 2 ;-axis parallel to the beam di¬ 
rection and X and y respectively are the horizontal and vertical coordinates 
perpendicular to the 2 ;-axis. The polar angle 9 is measured with respect to 
the z-axis and the azimuthal angle 0 is the angle with the x-axis in the xy 
plane. 

The GDG and VXDS give a combined momentum resolution of <Jp^/p± 
= 0.01 © 0.0026p±, where p± is the track momentum transverse to the beam 
axis in GeV/c. The VXDS consists of S barrels of Gharged Goupled Devices 
(GGD) at radii of 2.7, S.8 and 4.8 centimeters from the beam line, with 
S-hit solid angle coverage up to | cos6'|=0.85. The GGD pixels are cubic 
active volumes of 20 pm on each side, which provide precise SD spatial hits. 
This high granularity in SD space enable our track finding algorithm to limit 
the hit mis-assignment rate to only 0.2%. The spatial resolution on the hit 
cluster centroid achieved after a track based GGD alignment jHj is ~4 pm 
in both azimuthal and directions. The resultant tracking resolution for 
high momentum tracks, as measured from the miss distances of two tracks 
near the IP in ^ p+p“ events, is 7.7pm in r0 and 9.6pm in rz. The 
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multiple scattering contribution to the track impact parameter resolution 
can be approximately expressed as ^^.^ 3/2 g hm- These numbers are roughly 
a factor of two better than a typical vertex detector at LEP. This resolution 
advantage combined with the small and stable SLC IP information is crucial 
in establishing the feasibility of the measurement techniques, in particular 
for the case of the inclusive charm tagging for the Rc measurement. 

For the purpose of estimating the efficiencies and purities of the event 
selection and flavor tagging procedures, we use a detailed Monte Carlo (MC) 
simulation of the detector. The JETSET 7.4 [TU] event generator is used, with 
parameter values tuned to hadronic annihilation data nn, combined 
with a simulation of the SLD based on GEANT 3.21 P2]- The simulations 
of heavy hadron production and charm decays are described in [S]. The B 
decay simulation is adapted from the CLEO QQ MC with additional tuning 
by SLD (see appendix B of ini ) to match a wider range of T(4S') B decay 
data. 

For the hadronic event selection, we use a set of well-measured tracks, 
consistent with originating from within the beam pipe radius of 2.5 cm. The 
selected events have at least 7 tracks with p± >0.2GeV/c and within 5 cm 
from the interaction point along the beam axis. There are at least 3 tracks 
with two or more VXD hits. The event visible energy Eyis, calculated from 
charged tracks with the charged pion mass assigned, must be at least 18 GeV. 
To ensure the events are well contained within the detector hducial volume, 
we require the thrust-axis polar angle w.r.t. the beamline, 0 ^, calculated 
using calorimeter clusters, to be within | cos 6 *r| <0.7. We only include events 
with <3 jets for the analysis, where the jets are reconstructed from tracks 
using the JADE algorithm ^3] with i/cut = 0.02. This last requirement is 
to reject higher jet multiplicity events where the division of events into two 
hemispheres is no longer reliable for partitioning the events with one heavy 
hadron in each hemisphere. 

The selected event sample contains 191,770 events from the 1997-98 runs 
and 29,996 events from the 1996 run. We analyze the two data samples 
separately because the small 1996 sample, recorded during the early VXD3 
commissioning period, has a lower overall VXD hit efficiency due to electron¬ 
ics problems and some radiation damage. This effect was not present in the 
97-98 period when the VXD3 ran at ~10 degrees colder and the operation 
of the electronics was much more stable. 

The estimated background in the hadronic event sample is negligible. 


6 



The Monte Carlo hadrons event statistics used in the analysis have 

MC:data event ratios of 4:1, 15:1, and 22:1 for light flavor quarks {uds), cc, 
and bb events, respectively. 

3 Heavy Flavor Tagging 

Decays of bosons into charm and bottom quarks can be distinguished from 
those into light flavors {u,d and s) by searching for the heavy hadron decay 
vertices displaced from the event interaction point (IP). Because they are 
produced with high energy and have long lifetimes, heavy hadrons generally 
travel distances of millimeters before decaying. In a jet from a light quark 
the tracks will appear to come from one point in space, the event IP. In 
a charm jet some of the tracks may not point back to the IP, and if the 
charmed hadron decays into more than one charged particle there will be a 
secondary vertex (SV) in addition to the IP. Bottom jets will also exhibit 
secondary vertices, and if there are sufficient particles produced at the b and 
c decay points it is possible to hnd more than one displaced vertex. The 
reconstructed secondary vertices and their associated kinematic properties 
serve as the primary basis for flavor tagging for both bb and cc events in this 
analysis. 

3.1 IP Reconstruction 

To search for displaced vertices the position of the IP must be precisely 
known. In each event the position of the IP in the plane transverse to the 
beam axis is determined by htting all tracks that are compatible with coming 
from the IP to a common vertex. Because the SLC luminous region is small 
and stable in the xy plane, sets of 30 time sequential hadronic events are 
averaged to obtain a more precise determination of the xy IP position (details 
in appendix A of H)- The xy IP resolution is measured from the impact 
parameter distributions of Z^ —> events, deconvolved from the track 

impact parameter resolution measured from the two track miss distance, to 
obtain an IP resolution of 3.2 ym. 

Because the SLC luminous region is larger in z (700 /rm), the z position of 
the IP must be found event-by-event. Tracks with VXD hits are extrapolated 
to their point of closest approach (POCA) in xy to the precisely determined 
transverse IP position. Tracks with impact parameters of more than 500 /im 
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or 3 times the error (a) from the IP are excluded, and the z position of the IP 
is taken from the median z at POCA of the remaining tracks. The resolution 
of this method is found from the Monte Carlo simulation to be 10/11/17 pm 
for uds/cc/bb events. 

3.2 Secondary Vertex Reconstruction 

Secondary vertices are found using a topological algorithm QSI. This method 
searches for space points of large track density in 3 dimensions. Each track 
is parameterized by a Gaussian probability density tube f{r) with a width 
equal to the uncertainty in the track position at its POCA to the IP, ro: 

- (xq+ ^ - {zq + y tan X) \j^ | 

The hrst term is a parabolic approximation to the track’s circular trajec¬ 
tory in the xy plane, where k is a function of the track’s charge and transverse 
momentum and of the SLD magnetic held. The second term represents the 
linear trajectory of the track in the rz plane, where A is the track dip angle 
from the vertical. The a parameters are the uncertainties in the track posi¬ 
tions after extrapolation to Pq for the two projections. The function fi{f) is 
formed for each track i under consideration and used to construct the vertex 
probability function V (r). Also included is /o(T), a 7 x 7 x 20 pm {x x y x z) 
Gaussian ellipsoid centered at the IP position. 

Secondary vertices are found by searching for local maxima in V (r) that 
are well-separated from the peak at the IP position. The tracks whose den¬ 
sity functions contribute to such a local maximum are then identihed as 
originating from a secondary vertex (SV). 

A loose set of cuts is applied to tracks used for secondary vertex recon¬ 
struction. Tracks are required to have > 2 VXD hits and p± > 250MeV/c. 
Tracks with 3D impact parameter > 3 mm or consistent with originating 
from a or A decay, or 7 conversion are also removed. Each event is 
divided into two hemispheres using the thrust axis, and the vertexing pro¬ 
cedure is performed in each using only the tracks in that hemisphere. The 








identified vertices are required to be within 2.3 cm of the center of the beam 
pipe to remove false vertices from interactions with the detector material. A 
cut on the secondary vertex invariant mass M of \M — M^o | < 0.015 GeV/c^ 
removes any decays that survived the track cuts. The remaining vertices 
are then passed through a neural network uni to improve the background 
rejection further. The input variables are the flight distance from the IP to 
the vertex (H), that distance normalized by its error (D/an), and the angle 
(ppD between the flight direction D and the total momentum vector of the 
vertex P. These quantities are shown in Figure ^ along with the output 
value of the neural network {yvtx)- A good vertex in simulation is dehned 
as one which contains only tracks from heavy hadron decays, with no tracks 
originating from the IP, strange particle decays, or other sources. Vertices 
with y^tx >0.7 are retained. At least one secondary vertex passing this cut 
is found in 72.7% of bottom, 28.2% of charm, and 0.41% of light quark event 
hemispheres in the Monte Carlo. Around 16% of the hemispheres in b events 
have more than one selected secondary vertex. 

3.3 Track Attachment 

Due to the cascade nature of b hadron decays, tracks from the heavy hadron 
may not all originate from the same space point. Therefore, a process of 
attaching tracks to the secondary vertex (SV) has been developed to recover 
this information using a second neural network. The hrst four inputs are 
dehned at the point of closest approach of the track to the axis joining the 
secondary vertex to the IP. They are the transverse distance from the track 
to that axis (T), the distance from the IP along that axis to the POCA (L), 
that distance divided by the hight distance of the SV from the IP {L/D), 
and the angle of the track to the IP-SV axis (a). The last input is the 3D 
impact parameter of the track to the IP normalized by its error (b/at,). These 
quantities are shown schematically in Figure |21 The distributions are shown 
in Figure 121 along with the neural network output value {ytrk)- The network 
is trained to accept only tracks which come from a 6 or c hadron decay, 
and to reject tracks from the IP or from strange particle decays or detector 
interaction products. If more than one secondary vertex was found in the 
hemisphere the attachment procedure is tried for each track-SV combination. 
Tracks with ytrk > 0.6 are added to the list of secondary vertex tracks. This 
value is chosen to minimize the number of fake tracks being attached to 
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Figure 1: Distributions of seed vertex selection variables: (a) distance from IP 
D, (b) normalized distance from IP (c) angle between flight direction 

and vertex momentum (l)pD, (d) neural network output y^tx- A good vertex 
contains only heavy hadron decay tracks. The arrow indicates the accepted 
region. 
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Figure 2: Schematic illustration of the quantities used in the track- 
attachment procedure described in the text (not to scale). 



charm vertices, which than could mimic a b decay. 

3.4 Flavor Discrimination 

At this point, for each hemisphere there is a list of selected tracks. For 
hemispheres with no selected secondary vertices the list is empty, otherwise 
it includes the tracks in the secondary vertices and any cascade tracks which 
have been attached. From this list several signatures can be computed to 
discriminate between bottom/charm/light event hemispheres. These are the 
invariant mass of the selected tracks corrected for missing P* (M/jem), the 
total momentum sum of the selected tracks (Phem), the distance from the IP 
to the vertex obtained by htting all of the selected tracks and the 

total number of selected tracks {Nhem)- 

The four signatures given above are used as inputs for a neural network 
trained to distinguish hemispheres in bottom/charm/light events. The four 
inputs and the neural network output yhem are shown in Figure El The Pt 
corrected mass M/^em j2] , is a particularly powerful discriminator to separate 
bottom and charm. For a detector with high precision vertexing capability, 
decay vertices from b and c hadrons can both be very distinctively separated 
from the IP at high efficiency. The large b hadron mass is then the key 
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Tracks Tracks Tracks 


Figure 3: Distributions of the cascade track selection variables described in 
the text: (a) T, (b) L, (c) L/D, (d) a, (e) h/ab, (f) neural network output 
Utrk- A good track is one which originates from a heavy hadron decay. The 
arrow indicates the accepted region. 
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kinematic information to allow separation of bb and cc with high purity, 
which is crucial for these precision measurements. The raw vertex mass 
from the selected charged tracks can already give high purity b tags once 
requiring mass above the charm threshold, but many b decays with missing 
neutrals have an apparent low vertex mass. With the secondary vertex and 
PV positions very precisely measured at SLD, the b hadron flight direction 
can be derived to compare with the vector sum of the selected secondary 
track momenta to estimate the missing Pt with respect to the b hadron flight 
direction. A corrected vertex mass is then calculated to compensate for the 
derived minimal missing mass, taking into account the vertex position errors: 

Mfiem = \/ + |-Pt| (3) 

where rric is the invariant mass of the tracks attached to the secondary vertex 
in the hemisphere. 

Vertices in c quark jets, near the charm mass threshold typically have 
small missing P^, while many b vertices near the charm mass threshold receive 
a large missing Pt correction to become well separated from cc events. The 
missing particles for the low mass b vertices generally lead to a lower apparent 
momentum calculated from the visible charged tracks so that the correlation 
between Phem and Mhem presents another effective handle to further improve 
the flavor separation, as shown in (f) and (g) of Figure |3J comparing c and b 
jets. 

The flavor selection neural network is trained to put charm event hemi¬ 
spheres near yhem = 0, bottom event hemispheres near yhem = 1, and light- 
flavor background near yhem = 0.5. This allows a simple selection of charm 
(bottom) event hemispheres by specifying an upper (lower) limit for the out¬ 
put value yhem- Figure shows the ranges of purity vs. efficiency which can 
be obtained for charm and bottom event hemisphere tagging by adjusting 
this one cut. As can be seen in Figures [T] and 01 the inputs used for the 
neural networks at various stages of the yhem tag construction are in reason¬ 
able agreement between data and MC. However, an exact agreement is not 
essential as the main tagging efficiencies will be measured directly from the 
data. 
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Hemis Hemis Hemis 


Figure 4: Distributions of flavor discrimination variables: (a) Mhem, (b) 
Phem, (c) Dhem, (d) Nhem, (e) neural network output l/hem, (f) Phem vs. Mhem 
correlation for c jets (g) same for b jets. 
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Figure 5: Purity vs. efficiency for hemispheres in (a) charm and (b) bottom 
events, as the selection neural network cut is varied. 



b efficiency 



4 Measurement Method and Results 

Each event is divided into two hemispheres by the plane perpendicular to 
the thrust axis. Both measurements apply flavor tags on each hemisphere 
separately to derive Rf, and Rc, and measure the major tagging efficiencies 
from data simultaneously. 

4.1 Rb measurement 

For the Rb measurement, we dehne a single 6-tag with yhem >0.75 and apply 
the classical double tag method as in our previous measurement 0. By 
counting the fraction of tagged hemispheres, Fg, and the fraction of events 
with both hemispheres tagged. Fa, we can measure Rb by iteratively solving 
the following equations: 

Fg ^bRb ^cRc F Rb Rc) ^ 

Fd = Cb^^Rb F Cc^^RcF ~ Rh ~ Rc)- (4) 
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where eudsi^c and are the hemisphere tagging efficiencies for uds, c and b 

^double 

hemispheres, respectively, and C„ = —, are the hemisphere tag correla- 
tions for b and c events. We ignore the correlation for uds events since we 
expect 0.23 donble tagged events. A Standard Model valne of i?c=0.1723 is 
assnmed. The dependence on MC simulation is greatly reduced by measur¬ 
ing Cfe from data, while only the small ec,euds and correlations are estimated 
from MC. The measurement results and estimated MC parameters are tabu¬ 
lated in Tabled for the 96 and 97-98 data samples. The errors are statistical 
only for both measured and MC parameters. The Rb result has been cor¬ 
rected by —0.00149 for hadronic event selection bias and by -1-0.00024 for Zy 
interference effect. 



96 

97-98 

Rb 

0.21432 

± 

0.00289 

0.21624 

± 

0.00104 

efe(data)(%) 

56.63 

± 

0.67 

62.01 

± 

0.24 

efe(MC)(%) 

56.91 

± 

0.08 

61.78 

± 

0.03 

Hfe (%) 

97.8 

± 

0.4 

97.9 

± 

0.1 

ec(MC)(%) 

1.24 

± 

0.02 

1.19 

± 

0.01 

6„,,(MC)(%) 

0.113 

± 

0.006 

0.134 

± 

0.003 

Cfe-UMC) 

0.0067 

± 

0.0015 

-0.00012 

± 

0.00049 

a-l(MC) 

0.16 

± 

0.26 

0.30 

± 

0.12 

N events 

29996 



191770 



N single tag 

3315 



20738 



N double tag 

2091 



16048 




Table 1: Rb result and tagging performance parameters. life is the hemisphere 
6-tag purity. All errors are statistical only. 


4.2 Rc measurement 

Unlike the case of the Rb measurement where the single yhem > 0.75 tag has 
already collected b hemispheres with both high efficiency and high purity, 
the charm events are more spread out in the yhem distribution. Since the 
number of double tagged events is very sensitive to the tag efficiency, the 
lower efficiency for charm tagging makes the use of different tags with quite 
different purities more prohtable. We therefore use a multi-tag approach for 
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the Rc measurement. 

We divide hemispheres with a secondary vertex into four tag categories, 
depending on the output value of the neural network, i/hem- A b tag (tag 4) 
with i/hem > 0.75, a charm tag (tag 1) with i/hem < 0.3, a low-purity 6-tag 
(tag 3) with 0.5 < yhem < 0.75 and a low-purity charm tag (tag 2) with 
0.3 < yhem < 0.5. Hemispheres without a secondary vertex are classihed in 
the tag 0 category. In total we have 15 different event categories Nij, i < j < 
4 for the different tag combinations with a predicted fraction of the number 
of events Gij : 

+ ^AC:?Rc + -Rt- R,)] (5) 

with K = 1 for i = j and 2 for i ^ j. is the efficiency for quark q to give 
a tag i. is the tag correlation between tag i and tag j similarly dehned 
as in the Rb case. Since Nqq is determined by the total number of events 
and the other event counts, we have 14 independent equations for the event 
fractions Fij = A^p/A'totah 

A small number of uds events produce a secondary vertex. This has two 
reasons. The hrst is gluon splitting to bb and cc. The result is a real heavy 
quark of which the decay is well modeled in our simulation. These events 
mainly populate the high purity b and c tag categories and the errors on the 
measured rates dominate the systematic uncertainty assigned to this effect. 
The other source is incorrect reconstruction in our detector. These events 
typically cluster around yhem = 0.5. This is less well modeled, and in the 
next section we estimate the error on the rate predicted by our simulation to 
be 10%. The Nq 2 and A"o 3 categories contain most of these events. To avoid 
a bias from these low purity bins, we de-weight them by taking as the error 
on the event fraction, the binomial error on the bin contents and the 
systematic effect from varying the non-gluon splitting uds efficiency by 10%, 
added in quadrature. 

We minimize 

( 6 ) 

as a function of the following 9 parameters: R^ el (j = 1,4) and 

(j = 1,3). The c-quark efficiency for the 6-tag, e^, all light quark efficiencies, 
^iids (i = 1)4) and the hemisphere correlations C**! are taken from Monte 
Carlo. Only a few of the correlations are different from 1 in a statistically 
signihcant way. The others are set to 1. 
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The Rc fit results are summarized in Table |21 The Rc values have been 
corrected by —0.0004 for Z 7 interference, the event selection bias is zero in 
simulation. The Rc values are given at a central Rb value of 0.2157. The 
measured value of Rb agrees with the determination from the Rb measure¬ 
ment. There is a good agreement between the efficiencies in Monte Carlo 
and Data in the high purity tags 1 and 4. The efficiency for the lower pu¬ 
rity tags 2 and 3 for charm is higher in Monte Carlo than in the 97-98 data 
sample. The Monte Carlo efficiency is quite sensitive to physics parameters 
like the charmed hadron production fractions, their decay multiplicities and 
their lifetimes. The measured value of the efficiency can be reproduced in 
the Monte Carlo, by varying some of these parameters within their allowed 
range as is done in the study of systematic uncertainties. Since we extract 
the efficiencies from the data, the measured value of Rc is insensitive to these 
variations. 


5 Systematic Errors and Cross-checks 

The systematic uncertainties on Rb and Rc result from a combination of 
detector related effects and physics uncertainties in the simulation which 
affect our estimates of Cc, ends and Cg in the case of Rb and 64 , and all 
correlations for Rc- These systematic uncertainties are listed in Table [7| for 
the combined results for 96 and 97-98 data, which are analyzed separately 
initially. 

5.1 Hemisphere correlation cross check 

With the statistical precision down to the <0.5% level for Rb and the <2% 
level for Rc, the systematic uncertainties of the small subtle effects of hemi¬ 
sphere tag correlation become important. As the correlations have a variety 
of origins and the evaluation of their uncertainties will be spread over many 
subsections to follow, we will first discuss the correlation sources to estab¬ 
lish an understanding of the magnitude of their effects. Primarily as a cross 
check to constrain possible missing systematic sources for the hemisphere 
correlations, we decompose the efficiency correlation of the b and c tags into 
a set of approximately independent components which represent the major 
known sources of correlation between the two hemispheres in the bb and cc 
MC events. 
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96 

97-98 


data 

MC 

data 

MC 

Rc 

0.1678 ±0.0091 


0.1752 ±0.0033 


el (%) 

2.22 ±0.18 

2.01 

2.52 ±0.07 

2.36 

6 ^ (%) 

16.09 ±0.90 

15.75 

17.95 ±0.36 

18.56 

4.. (%) 


0.088 


0.063 

el (%) 

2.83 ±0.20 

2.56 

2.97 ±0.07 

3.00 

(%) 

4.93 ±0.36 

5.05 

5.02 ±0.13 

5.72 

elds (%) 


0.117 


0.124 

e! (%) 

4.68 ±0.24 

4.40 

5.10 ±0.09 

4.98 

6 ^ (%) 

2.46 ±0.27 

2.59 

2.23 ±0.10 

2.87 

elds (%) 


0.150 


0.140 

et (%) 

56.8 ±0.64 

56.9 

62.18 ±0.23 

61.78 

4 (%) 


1.23 


1.19 

4.. (%) 


0.113 


0.134 

- 1 


-0.0068 


0.012 

-1 


0.0032 


0.0011 

Cf-l 


0.0064 


-0.0015 

Cf - 1 


0.0067 


- 0.00012 

other CP — 1 
xVd.o.f. 

6.66/5 

0 

7.05/5 

0 


Table 2: Rc results and tagging performance parameters. The errors on the 
MC efficiencies are small compared to systematic uncertainties and therefore 
omitted 


To focus on understanding the physics sources, we use the large uni¬ 
form 97-98 MC sample without tracking efficiency and resolution corrections, 
to avoid the statistical fluctuations introduced by the tracking corrections. 
The overall tracking systematic effects will be evaluated separately in sec¬ 
tion EISUl For the Rc measurement, there are many correlations between the 
different tags. We present a representative c-tag of Hhem < 0.4, for similar 
set of effects as for the 6 -tag. For this study we need to identify the relevant 
kinematic and geometric variables to see how the tagging efficiencies depend 
on them and how the two hemispheres correlate on these variables. 

The primary vertex (PV) shared between the two hemispheres is an ob¬ 
vious source of correlation. A misreconstructed PV results in a negative 
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correlation if the displacement is along the thrust axis. Displacement of the 
PV transverse to the thrust axis would tend to positively correlate the two 
hemispheres. The small and stable SLD IP in the xy view is an average 
beam position over many events which greatly reduces the chance of large 
PV displacement in xy, and restricts the remaining effect to be only through 
the IP coordinate, which is reconstructed event by event. The total PV 
effect is studied by simply comparing results using the reconstructed PV and 
using the MC truth PV. 

Another major source is geometric correlations. The two hemispheres 
are exactly back to back by dehnition. Due the cylindrical geometry of 
the detector, larger | cos 6^1 typically means more multiple scattering, worse 
tracking efficiency and resolution, worse radial alignment etc. This will affect 
the hemisphere at cos 6 and the opposite hemisphere at — cos 6 in the same 
way. We calculate the tagging efficiency in bins of the hemisphere axis cos 6 
bin and construct a function e(cos6'). The cos6 correlation component can 
then be estimated from 

C = ^ e(cos 6i) ■ e(— cos 6i) /(cos 6i) 

i 

where i is the index of the cos 9 bin, / is the fraction of all events in bin i 
and e is the average tagging efficiency of all hemispheres. 

Although the cylindrical geometry of the detector ensures uniformity at 
first order in azimuth, local performance variation of detector elements can 
still introduce an effect. Similar to the 9 component analysis, the / compo¬ 
nent is estimated from 


^forward(^4^i) ' ^backward{,4^i T f /C 

i 

where the efficiency parameterization is done for forward {+z) and backward 
{—z) hemispheres separately. 

Tagging correlation due to variation of detector performance in time also 
needs to be addressed. This is in fact the primary reason that the 97-98 data 
and 96 data are analyzed separately. There is a signihcant VXDS perfor¬ 
mance variation during the 96 run, and the Monte Carlo simulation repro¬ 
duces these time dependent effects in detail. The time dependent component 
is estimated from 
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where Ri is a group of runs adjacent in time, f{Ri) is the fraction of all 
events in this group of runs. 

The h/ c-tag efficiencies have a signihcant dependence on the heavy hadron 
momentum {Ph) and its angle to the hemisphere axis (.^). There are various 
causes, such as gluon radiation, which can result in correlations of the heavy 
hadron momenta and angular distributions. For the momentum correlation 
component, we parameterize the hemisphere tagging efficiency in a 2D grid 
of e{PH,0 then sum over all events i for the correlation for hemispheres 1 
and 2: 

i 

This component is only evaluated for the normal cases where two heavy 
quarks are in opposite hemispheres. The joint effect of the heavy hadron 
momentum and thrust angle is ~20% larger than the effect from momentum 
alone. 

For the extreme case of hard gluon radiation, with two heavy quarks 
recoiling into the same hemisphere opposite the hard gluon, we simply calcu¬ 
late the effect if including and excluding these events in the overall sample. 
The magnitudes of the various components and their sum compared with the 
overall correlations calculated from the MC hemisphere and double tag rates 
are shown in Table El An important observation is that the magnitude of all 


Component 

{Cb-tag 

97-98 

- 1) X 10^ 
96 

(Cc-tag 

97-98 

- 1) X 10*^ 
96 

Primary vertex 

+46 

+13 

-480 

-640 

Geometrical correlation 6 

+49 

+60 

+120 

-120 

Geometrical correlation (j) 

-4 

+212 

-10 

+170 

Time dependence 

11 

+434 

+40 

+870 

B/D momentum and thrust angle 

+107 

+95 

+820 

+760 

Hard gluon radiation 

-37 

-23 

+510 

+410 

Component sum 

+170 

+670 

+1000 

+1450 

MC overall correlation 

+42 

+891 

+420 

+1510 

MC statistical error 

±47 

±113 

±390 

±1210 

discrepancy 

+ 128 

-121 

+580 

-60 


Table 3: Hemisphere correlation component check results for 6-tag and c-tag. 


components are small, which is especially true for the 6-tag. The change in 
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g ^ cc (%) g —^bb (%) 

LEPEWWG standard 
JETSET 

2.96T0.38 0.254T0.051 

1.357 0.142 


Table 4: The world average measurements of g —y QQ compared to the pre¬ 
dictions of the JETSET generator. 


the value of the correlation C translates to the fractional error of Rb and Re- 
So effects at typically 0.1% (1%) or less for the 6(c)-tag can be compared with 
the Rb (Re) statistical error of ~0.5% (2%). This illustrates the importance 
of the large tagging efficiency in driving a much smaller correlation and re¬ 
ducing sensitivity to correlation uncertainties, which is even apparent when 
comparing the 6-tag with the c-tag. In the limit of 100% tagging efficiency, 
the correlation becomes irrelevant. 

The 96 samples are not large enough to draw detailed conclusions given 
the large MC statistical error, but the difference of the overall correlation 
between 96 and 97-98 is well explained by the known time and 0 dependent 
effects in 96. For the more statistically significant test with the 97-98 sample, 
there is a noticeable discrepancy between the 6-tag component sum and the 
MC overall correlation. Although at least part of the discrepancy could be 
statistical for this approximate cross check, we take half of discrepancies from 
the 97-98 analysis for both the 6-tag and c-tag as an additional systematic 
error associated with possible missing sources. 

5.2 Physics systematics 

The physics systematic errors are mostly assigned by reweighting the nominal 
simulated distributions to an alternative set of distributions which correspond 
to the world average measurements and uncertainties of the underlying MC 
physics parameters in|. 

5.2.1 g QQ effect 

The g ^ bb and g ^ cc production rates are varied according to the experi¬ 
mental averages in|. These measurement rates are significantly higher than 
the default JETSET MC as shown in Table 01 The main sensitivity to the 
gluon splitting uncertainties is through the uds tagging efficiencies. Its ef- 
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Figure 6: Components of the Rf, systematic error as a function of the flavor 
separation yhem cut for the various g —>• QQ processes. 


feet on charm background under the &-tag is also noticeable. In the case of 
the Rh measurement, we can see in detail in Fig. IHlhow g —>■ QQ affects the 
systematic uncertainty on as a function of yhem cut. 

5.2.2 uds physics modeling 

Among the fake uds tags, a significant fraction involve and A decay 
product. The MC generator level and A production are varied by ±10%, 
as recommended in HZI, to derive the systematic uncertainty due to light 
hadron production. 
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5.2.3 charm production and decay modeling 

The charm production and decay modeling affect the Rb and Rc measure¬ 
ments through the uncertainty on the small charm tagging efficiency in the 
6 -tag. They also affect the Rc measurement through the charm tag hemi¬ 
sphere correlations. Many of the production and decay uncertainties only 
affect charm tag correlations indirectly, e.g. through the primary vertex re¬ 
construction. 

The different charmed hadrons have very different lifetimes and also 
rather different decay multiplicities so that the tagging efficiencies are also 
rather different. The D^, D^, Dg and lifetimes are varied according to the 
errors on their measured values as in PDG QSI. The production fractions of 
these different charmed hadrons are varied according to the recommendations 

of im. 

The tagging efficiency is in general higher for higher momentum charmed 
hadrons so that the energy spectrum of the charmed hadrons is another 
source of uncertainty. Charmed hadron fragmentation is studied by varying 
the average scaled energy in the range {xe) = 0.484T0.008 using the Peterson 
fragmentation fnnction as well as by comparing the difference between 
the Peterson and Bowler models [20] for the same values of {xe)- 

Charm tagging efficiency is sensitive to the charged decay mnltiplicity. 
The fraction of decays with fewer than two tracks is a crucial source of in¬ 
efficiencies. Higher multiplicity vertices are easier to identify, but the softer 
decay product momenta do not necessarily make the vertex resolution better. 
The uncertainty due to D decay charged multiplicities is estimated by vary¬ 
ing each decay multiplicity fraction according to the Mark-IIl measurement 
errors EH. for and Dg in turn, with a specific scheme as described 

in jT7j. 

The production of in charm decays is another relevant source of uncer¬ 
tainty, although it may be somewhat correlated to the charged multiplicity. 
In the case of or all neutral decays, it is a significant source of charged 
multiplicity loss. In the case of tt+tt” decays, the decay product may 

affect the vertexing of charm decay differently depending on the decay 
length. There is no explicit recommended scheme from LEPEWWC. We 
reweight the multiplicity associated with charmed hadron decay to cor¬ 
respond to the average multiplicity range measured by Mark-III pT] . 

Charmed hadron decays with fewer neutral particles have higher charged 
mass and are therefore more likely to be mistagged as a b. Thus, an additional 
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systematic uncertainty is estimated by varying the rates of charmed hadron 
decays with no tt^s by ±10%. This is an SLD specihc estimate with particular 
relevance to mass tags, which is typically not included in LEP measurements. 

5.2.4 b production and decay modeling 

The b hadron production and decay modeling uncertainties only enter via the 
6 -tag hemisphere correlation. Since the B decays, and to a large extent the 
b jet fragmentations proceed independently between the two hemispheres, 
most of these effects (apart from the QCD gluon radiation effect discussed 
separately in the next section) enter indirectly through such effect as primary 
vertex reconstruction. 

The 6-hadron lifetimes are varied by a typical current measurement er¬ 
ror of ±0.05ps. We also vary the A;, production ratio (not included in the 
LEPEWWG recommendation) by current measurement uncertainty to ac¬ 
knowledge the fact that the A;, lifetime is signihcantly shorter than that of 
the B mesons. 6 hadron fragmentation is studied by varying the average 
scaled energy in the range {xe) = 0.702 ± 0.008 using the Peterson frag¬ 
mentation function as well as by comparing the difference between the 
Peterson and Bowler models ra for the same values of {xe)- The B decay 
charged multiplicity distributions are reweighted to correspond to an average 
charged multiplicity uncertainty of ±0.35. 

5.2.5 Gluon radiation effects on tag correlations 

Gluon radiation is a more direct source of correlation as it tends to simulta¬ 
neously lower the momenta of both heavy quarks as well as changing their 
directions from the back to back topology. The tags are generally sensitive 
to heavy hadron momentum and their directions w.r.t. the thrust axis. 

We hrst discuss the extreme case of a very hard gluon radiation which 
causes two heavy hadrons to recoil into the same hemisphere. This creates 
an anti-correlation of tagging efficiencies between the two hemispheres. Hard 
gluon radiation resulting in two H’s in the same hemisphere occurs at a rate 
of 2.45% in our simulation. This effect is reduced to 1.62% in the 66 MC 
passing the analysis event selection. In the case of the 6-tag, the two H’s 
have lower B momentum and wider angles to the thrust axis, so that the 
6 -tag efficiency is only slightly higher than a normal single 6 hemisphere. 
The c-tag is more sensitive to the reduced c-hadron momentum and the 
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more confusing kinematic situation results in a larger effect to the extent 
that the hemisphere with the two c-hadron even has a signihcantly lower 
efficiency than normal. For the actual systematic uncertainty, we follow the 
LEPEWWG recommendation to vary the MC rate by ±30%. A cross check 
measurement (Appendix B.4.1 in of the rate of hard gluon radiation 
is performed with 6-tags applied to each jet in 3 jet events. This analysis 
measured the ratio of same-hemisphere bb rates between data and MC to be 
0.82±0.09, well within the ±30% variation. When the tagging efficiency of 
the hemispheres with two S’s is very close to the normal hemispheres with 
one 5, there is a compensating effect between the double tag and hemisphere 
tag such that the overall correlation becomes very insensitive to the fraction 
of hard gluon radiation events. This results in a rather small error of SRb = 
—0.00002. In the case of Rc analysis, the hard gluon events in both bb and 
cc are weighted up by 30% for all tags and the combined effect of 6Rc = 
±0.00026 is dominated by hard gluons in cc events. 

A systematic uncertainty is also assigned to the momentum correlation 
between the two heavy hadron hemispheres, mostly due to gluon radiation 
and fragmentation effects, which in turn translate to a tagging efficiency cor¬ 
relation. In the case of the B momentum correlation in bb events, this is 
estimated by comparing the B momentum correlation in the HERWIG | 22 ] 
and JETSET jTUj event generators. At the parton level, all generators give 
a similar correlation of ~1.4% between the b quarks. At hadron level, the 
correlation coefficient for the 6-hadron momenta in the two hemispheres is 
1.55% in JETSET, while the largest deviation among different models is seen 
in HERWIG which gives up to ±0.8% higher B momentum correlation. We 
use half of this difference of 0.4% as the variation to estimate the system¬ 
atic error, according to the recommendation in nn. Given that different 
event selections can change the absolute correlation (excluding case of two 
Bs in same hemisphere, JETSET correlation reduces to 1.23%, and applying 
Njet < 4 cut further reduces this to 0.85%), the ratio of 0.4%/1.55%=0.26 is 
taken as the fractional uncertainty on the 6 momentum correlation effect. A 
cross check analysis (Appendix B.4.2 in of B hadron energy correlation 
is performed using the observed vertex momentum from charged tracks in 
double tagged events. The B momentum correlation measured is verihed to 
be in good agreement with MG to ±20%. 

The angular distribution of B flight direction and thrust axis is also 
checked to be in good agreement between data and MG, where the B direc¬ 
tion is approximated by the line joining 6-tag vertex and PV. The component 
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of the b tagging correlation due to the B momentum correlation is estimated 
to correspond to Cf, — 1 = 0.00107, as is described in section HTTl which trans¬ 
lates to an error on Ri, of 0.00006. There is no equivalent recommendation in 
UTI for c-hadron momentum correlation in cc events, mainly because this is 
only relevant for double tag analysis which is not done at LEP. We similarly 
take ±26% of the H-momentum correlation component in cc correlation as 
a systematic uncertainty, which corresponds to 6Cc-tag = 0.0020 and 6Rc of 
0.00022. 

5.3 Detector systematics 

5.3.1 Tracking resolution and efficiency 

The tracking resolution systematic effects, primarily due to residual detector 
misalignment are estimated from the observed shifts in the track impact 
parameter distributions as a function of 0 and 6 in both r-0 and r-z planes. 
The typical impact parameter biases observed in the data are ~ 2.5 pm 
(~ 5 pm) in the r-0 (r-z) plane. A correction procedure is applied so that the 
MC tracks match the mean bias values of the data in various (0, 6) regions. 
This is a more realistic evaluation of alignment bias effects, where tracks 
passing the same detector region are biased in a correlated manner. The 
actual corrections are implemented by dividing 0 into 40 regions according 
to VXDS ladder triplet boundaries and 4 sections in 6. A detailed description 
of the resolution corrections can be found in Half of this correction is 
taken as the variation to evaluate the track resolution systematic uncertainty. 

The uncertainty in the tracking efficiency is evaluated from a comparison 
between data and MC for the fraction of all CDC tracks which pass a set 
of quality cuts, and the fraction of good CDC tracks extrapolating close 
to the IP that does not have associated VXD hits. These studies indicate 
that the MC overestimates the tracking efficiency by ~1.5% on average. A 
procedure for the random removal of tracks in bins of pt, 0, and 6 is used 
to correct the MC for this difference. Secondary vertex charge distributions 
are used as an independent check to verify that they agree better between 
data and MC with these corrections applied, as shown in Fig. [Tj This tests 
directly the relevant secondary vertex tracks independent of the modeling of 
fragmentation track multiplicity. The |Q| > 1 bins are particularly sensitive 
to tracking inefficiencies, while only weakly depend on heavy hadron decay 
multiplicity and production fractions. 
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Figure 7: Top: vertex charge distributions for yhem > 0.75 &-tag (left) and 
Uhem < 0.4 c-tag (right) hemispheres. The data (dots) are compared with 
MC predictions with or withont tracking corrections. Bottom: the ratio of 
MC over data for before and after corrections 













5.3.2 IP tail 


Another potential systematic source is the possible tails in the IP position 
determination which can cause large tagging asymmetries. The ii~ 

events have no indication of systematic tail effects. A further study [22] is 
done by examining the distribution of the distance between the hemisphere 
axis drawn through the htted event primary vertex in r — 0 to the IP for 
2-jet hadronic events with no secondary decay tracks. 

There is no significant discrepancy between data and MC for 97-98 data, 
but there is a detectable tail in the 96 data with ~100 /rm width for 0.5% of 
the events. The full effect of the tail is treated as a systematic uncertainty 
for the 96 data. 

5.3.3 Geometric and temporal correlation effects 

Besides the issue of the uncertainty in the overall tracking efficiency, MC 
modeling uncertainties on the efficiency nonuniformity due to local detector 
inefficiencies and time variations can lead to additional uncertainties on tag¬ 
ging correlations. The contribution of the 9 and (j) tagging efficiency depen¬ 
dence to the overall tagging efficiency correlation are estimated in sect ion UnTl 
We estimate the uncertainty of the MC modeling of 6 and 0 dependency, by 
comparing the hemisphere tagging fractions between data and MC for the 
yhem > 0.75 6-tag and representative yhem < 0.4 c-tag. There is a gener¬ 
ally good agreement between the data and MC within the statistical errors. 
The only significant effect, the variation in 0 for 96 data, due to electronics 
failures for some of the periods, is reasonably simulated by the MC. To es¬ 
timate the correlation uncertainty, we simply take the ratio of the all input 
hemisphere tagging fractions between data and MC for each 0 and 6 bin, to 
reweight the MC true signal tagging efficiency for the corresponding bin, as 
an approximated deviation. The reweighted 0 and 6 dependent efficiencies 
are then used to recalculate the tagging efficiency correlation components 
as described in section o The change of the tagging correlation from the 
reweighting is taken as a systematic error. Note that correlation component 
calculations are insensitive to overall efficiency differences between data and 
MC which may enter in the reweighting, but only to shape variations. The 
resulting systematic errors on the geometrical correlations are summarized 
in Table El The correlation systematic error on ACb-tag translates directly 
to the fractional error 6Rb/Rb, and the error on the representative c-tag of 
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97-98 96 

97-98 96 

6 ehect 

0 ehect 

time dependence 

0.00016 0.00181 
0.00027 0.00037 

- 0.00029 

0.00043 0.00337 
0.00030 0.00409 

- 0.00055 


Table 5: Correlation systematic uncertainties due to geometrical and time 
dependent effects. 


ACc-tag translates approximately to the fractional error 6Rc/Rc- The use 
of data vs MC ratio unfortunately introduces statistical fluctuations on the 
estimated correlation change. A toy MC study indicates for example that 
the statistical fluctuation expected for the 0 effect estimate is ~0.00009 for 
97-98 6-tag and as large as ~0.0031 for the 96 c-tag. Some of the changes are 
consistent with statistical fluctuation, but we still conservatively take them 
as systematic uncertainties. 

For the time dependent effects on the tagging correlations there is a very 
good uniformity in the 97-98 data. The smaller 96 data sample does have 
signihcant time dependent effects, which contribute to the overall uncertain¬ 
ties. The resulting systematic errors on the temporal correlations for the 96 
analysis are summarized in Table El 

5.4 Event selection bias 

The hadronic event selection flavor bias is evaluated from the MC for the 
basic hadronic event selection procedure and the last step of Nj^t <4 cut 
separately. The event selection efficiencies and resulting bias on Rb and 
Rc are tabulated in Table El The basic hadronic selection passes a slightly 
higher fraction of bb events, which is expected from the known higher charged 
multiplicity and other observed kinematic differences of 6-jets compared to 
uds. Given that the effect is only a few times the statistical error a, we 
simply take the MC statistical error as an uncertainty for this stage. The 
Njet < 4 cut on the other hand has a more statistically signihcant effect. 
The rate of cc events passing the Nj^t < 4 cut is consistent with that for uds 
(most of the ARc is actually the compensating effect of Rb bias). 

The effect of the 6 quark mass on the > 4 jet rate has signihcant theo¬ 
retical uncertainties. The JETSET 7.4 MC used for our analysis is in fact 
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Stage 

Basic Had. Sel. Nj^t < 4 

Total 

uds efficiency 
c efficiency 
b efficiency 

58.330 ±0.050% 91.760 ± 0.037% 

58.592 ± 0.045% 91.695 ± 0.033% 

58.694 ± 0.036% 92.376 ± 0.025% 

53.524 ± 0.051% 
53.725 ± 0.045% 
54.220 ± 0.036% 

ARb 

ar. 

0.00089 ± 0.00016 0.00117 ± 0.00007 

0.00040 ± 0.00018 -0.00035 ± 0.00008 

0.00207 ±0.00017 

0.00005 ±0.00020 


Table 6: Event selection efficiencies and flavor bias for the uncorrected MC. 


known to have excessive suppression of gluon radiation for heavy quarks 
PK] . when comparing with data on 3 jet rates as used in the running b quark 
mass measurements [21]. We use the DELPHI measurement [211 of the ratio 
= fi{hb)/f\{uds) to evaluate the event selection bias in our analysis for 
the < 4 jet requirement, where denotes the fraction of events with > 4 
jets. The DELPHI measurement using the Cambridge jet finder on all final 
state hadrons gives R’^ = 0.89 ± 0.02 at j/t.ut=0.006, with only very small 
dependence. The overall 4 jet rate at this is very similar to the 4 jet rates 
using JADE Yclus jet finder on charged tracks with y^^t=0.02, as in our anal¬ 
ysis. We use the PYTHIA 6.228 generator 1271 as an intermediate reference 
to compute correction factors based on the ratios of (1 — R^) with the same 
jet finding algorithms. We obtain a scaling factor of S' = 0.50 ± 0.13 to be 
applied to the raw 4-jet cut bias in the JETSET 7.4 MC. The uncertainty 
consists of equal contributions from the DELPHI measurement uncertainty 
and from the observed difference between using charged tracks and all final- 
state hadrons as input to the jet finder algorithms. Applying this scaling 
factor to the raw 4-jet bias of ARb = 0.00117 in Table jH we obtain a cor¬ 
rected 4-jet bias of ARb = 0.00059 ± 0.00015. The event selection bias for cc 
is consistent with zero and therefore no correction is applied for R^. 

As a crosscheck to verify the Rb event selection bias estimate for the 4- 
jet cut, we also performed the analysis without the Nj^t < 4 cut for the 
larger 97-98 data sample. Taking into account the additional sensitivity to 
g —>• QQ, we And the difference in measured Rb due to the removal of < 4-jet 
requirement is 

Rb{no < 4 jet cut) — i?b(nominal analysis) = -1-0.00042 ± 0.00030 ± 0.00015 

where the first error is the uncorrelated statistical uncertainty and the second 
error is the systematic error on the 4-jet cut bias as estimated above. The 
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difference is consistent with zero within statistics. 


5.5 Result stability checks 

One possible source of systematic uncertainty is poor modeling of the uds 
background that gives a secondary vertex due to badly reconstructed tracks. 
As a cross-check we tried to extract this effect from the data. We dehne a 
uds tag by requiring no secondary vertex and no track with a normalized 3D 
impact parameter of more than 2a. This tag identihes about 50% of the uds, 
15% of the c and 1.4% of the b hemispheres. Adding this tag to the other 
ones, we £t for the same efficiencies as before plus the ratio of uds giving a 
secondary vertex in data over Monte Carlo. We find a value of 1.1 ± 0.1 for 
this ratio. 

We also examine the variation of the Rb result as a function of the mini¬ 
mum yhem cut for the b tag, over a wide range of yhem cut, as shown in Fig.|Hl 
The Rb result is stable within the la total uncertainty envelope. 

Since the Rc measurement uses the full range of yhem values, there is no 
equivalent plot for this measurement. We do check the stability of the result 
by moving the boundary between the high purity charm tag (1) and the low 
purity charm tag (2) from 0.15 to 0.45 in increments of 0.05. We find that the 
Rc result is stable within 0.0005. The low purity tags showed some difference 
in efficiency for charm events between data and Monte Carlo. Although we 
can reproduce the measured efficiency by varying some of the MC parameters 
in their allowed range we perform an additional cross-check by removing the 
low-purity tag 3 from the £t and re-measuring Rc. We hnd an insignihcant 
shift of +0.0009 ± 0.0014 in 97-98 and -0.0022 ± 0.0032 in 96. 

6 Conclusions 

We have measured the hadronic branching ratio of the to b quark and 
c quark with our 96-98 dataset of 400,000 hadronic decays. The 96 
and 97-98 results are combined, with all common systematic uncertainties, 
including detector uncertainties, treated as fully correlated. The combined 
results are: 

Rb = 0.21604 ±0.00098(stat) ±0.00073(sys) +0.00012(72^) 

Rc = 0.1744 ±0.0031(stat) ±0.0020(sys) +0.0006(72f,) 
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Figure 8: a) Tagging efficiencies and b purity vs yhem cut. b) variation of 
the measured value of Rb, with the inner error bar being the uncorrelated 
statistical error w.r.t. the nominal cut and the outer bar total error, c) 
Rb statistical and systematic error components as a function of the 6-tag 
selection cut, for the 97-8 analysis. The net uncertainty does not include 
event selection uncertainties. qq 













































For the Ri, measurement, combining this new measurement with our previ¬ 
ously published result on the 93-95 data j2], we obtain: 

Rb = 0.21594 ± 0.00094(stat) ± 0.00074(sys) ^ 0.00012(7?^) 

The relative weights in the combined average for the 93-95:96:97-98 (96:97- 
98) Rb {Rc) measurements are 7:9:84 (11:89), dominated by the 97-98 result. 
These measurements are in good agreement with the Standard Model ex¬ 
pectation of i?ft=0.2156 (for mi=178 GeV/c^) and i?c=0.1723. They can be 
compared with the average of LEP measurements [2H] from a total of ~16 M 
hadronic decays: 

Rb = 0.21643 ± 0.00073 
Rc = 0.1691 ± 0.0047 

In conclusion, we have exploited the high resolution vertexing capability and 
the small and the stable SLC IP for a precision test of Standard Model 
through the measurements of heavy quark production fractions in decays. 
Our new Rc result is by itself more precise than the current world average 
|E|. These measurements conhrm the Standard Model predictions at ~0.6% 
precision for Rb and 2.1% precision for Rc. 
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of Energy and National Science Foundation, the UK Particle Physics and As¬ 
tronomy Research Gouncil, the Istituto Nazionale di Fisica Nucleare of Italy 
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5Rb{lO-^) 

<5i?c(10-°) 

MC statistics 

13 

91 

g^hh 0.254±0.051% 

-24 

9 

g^cc 2.96±0.38% 

-23 

-101 

long lived light hadron production ±10% 

-1 

-1 

production 0.233 ± 0.028 

-10 

-6 

Ds production 0.102 ± 0.037 

-11 

-15 

c-baryon production 0.065 ± 0.029 

-11 

22 

charm fragmentation 

-18 

18 

lifetime 0.415 ± 0.004 ps 

-3 

8 

lifetime 1.057 ± 0.015 ps 

-2 

5 

Dg lifetime 0.467 ± 0.017 ps 

-3 

-3 

Ac lifetime 0.206 ± 0.012 ps 

-1 

-91 

D decay multiplicity 

-27 

60 

D decay 77° 

19 

56 

D decay no-7r° 

-9 

12 

B lifetime ±0.05 ps 

0 

5 

B decay (Nch) = 5.73 ± 0.35 

-20 

3 

b fragmentation 

4 

26 

Afe production fraction 0.074 ± 0.030 

5 

-2 

QCD hemisphere correlation 

6 

22 

hard gluon radiation 

-2 

26 

tag geometry dependency 

9 

17 

tag time dependency 

1 

1 

component correlation 

14 

45 

tracking resolution 

27 

22 

tracking efficiency 

13 

3 

{IP)xy tail 

2 

0 

event selection bias 

17 

20 

4 jet rate in b events 

15 

0 

7?c = 0.1723 ±0.0037 

Rb = 0.2157 ±0.0010 

-12 

-62 

Total (excl. Rb/c) 

73 

200 


Table 7; Summary of systematic uncertainties for Ri, and Rc- A — sign for 
an error means the value for Rq goes down when this parameter is varied 
upward. The errors are assumed to be symmetric. 
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